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ABSTRACT 



We report SOFIA/GREAT, HerschelfHWI, and ground-based velocity-resolved spectroscopy of carbon monoxide (CO) rotational 
transitions from J=2-l to J=16-15 toward two positions in the circum-nuclear disk (CND) in our Galactic center. Radiative transfer 
models were used to derive information on the physical state of the gas traced by CO. The excitation of the CO gas cannot be explained 
by a single physical component, but is clearly the superposition of various warm gas phases. In a two-component approach, our large 
velocity gradient (LVG) analysis suggests high temperatures of ~200 K with moderate gas densities of only ~10 45 cirr 3 for the bulk 
of the material. A higher excited phase, carrying ~20-30% of the column densities, is warmer (~300-500 K) but only slightly denser 
(~10 5 3 cm -3 ). These densities are too low to self-stabilize the clumps against their high internal turbulence and fall below the Roche 
density (>10 7 cm" 3 ) at 1.5 pc galactocentric distance. We conclude that the bulk of the material in the CND is not organized by 
self-gravity nor stable against tidal disruption, and must be transient. 
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1. Introduction 

At a distance of only 8.5 kpc, our Galactic center (GC) is a 
unique laboratory for studying the physical processes that also 
occur within extragalactic nuclei in general. The innermost cen- 
tral parsecs contain all the ingredients and processes that also af- 
fect the gas in other nuclei that are too distant for detailed study. 
Here we are most concerned with the massive gas clouds that lie 
within a few parsecs of the Galactic center's black hole, SgrA*. 
It is only in our own Galaxy's center that we can spatially resolve 
and study these likely reservoirs of matter that fuel the center's 
episodic activity. 

The global distribution of the interstellar matter in the very 
center of the Galactic nucleus is a — 1-1.5 pc radius ionized cav- 
ity centered on SgrA*. The cavity's outer boundary forms the 
sharp inner edge of the circum-nuclear disk (CND), a thin disk- 
like structure composed of dense molecular gas in filaments and 
streamers that extend to radius ~5 pc. Ionized material is present 
also in the cavity in the form of streamers with a "mini-spiral" 
shape. These streamers have high radial velocities and carry ma- 
terial deep into the central cavity and close to SgrA* (for more 
detailed references we refer to the review by Genzel et al. 20 10). 

Despite many investigations, the physical characteristics of 
the neutral gas remains a topic of debate. Estimates of the CND's 
molecular gas mass range between 10 4 (from dust) and 10 5 ~ 6 M 
(derived from gas tracers). Depending on the probe, density esti- 
mates for the dense clumps are 10 5 ~ 8 ctrT 3 and gas temperatures 
run from 50 to a few hundred K. The range of physical param- 
eters leaves open many questions about the nature and fate of 
the CND, for instance whether the clumps' densities exceed the 
Roche limit and make them stable against tidal disruption, or 
whether they are fully transient features. 



A prerequisite for a quantitative comparison between these 
competing scenarios is an in-depth investigation of the temper- 
ature and density distribution in this clumpy, filamentary ring. 
Early work from Genzel et al. (1985) detecting CO(16-15) with 
the KAO, cle arly indicated the ex istence of a gas phase with high 
temperatures. Harri s et al. 1(119 85) tied low-, mid-, and high-J ob- 
servations together for a full CO ex citation analysi s that gave 
Tktn -300 K and n(H 2 ) ~10 45 cm" 3 . iBradford et al.1 (120051) . re- 
visiting these data with new CO(7-6) observ ations derived warm 
200-300 K and moderate densities 5-7xl0 4 . IWhite et al. I d2003h 
used infrared ISO/LWS data to find gas with T kin -900 - 1400 K 
andn(H 2 )~ 10 4 - 10 4 " 



k5 cm -3 



In this contribution we present initial results of a system- 
atic and comprehensive investigation of the CND, combining 
velocity-resolved heterodyne spectroscopic observations with 
the far-infrared spectrometers GREAlQ on board of SOFIA and 
HIFI/PACS on Herschel, supplemented by data of lower-J CO 
transitions from APEX and the IRAM-30m telescope. 

Here we perform an in-depth CO excitation analysis toward 
the prominent northern and southern lobes of the CND (Fig. 1). 
The CO spectral energy distribution is from line-resolved obser- 
vations of most transitions of the CO rotational ladder up to J = 
16-15, and additional 13 CO isotopologues to confine the optical 
depth of the CO lines. 



1 GREAT, the German REceiver for Astronomy at Terahertz fre- 
quencies, is a development by the MPI fiir Radioastronomie and 
the KOSMA / Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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Table 1. Observed integrated line intensities 



R.A. offset [orcsec] 



Fig. 1. CND in the emission of CO(6-5), observed with the 
CHAMP + receiver at APEX. The CO intensity (color map in 
K km/s, in T m b) was integrated between -150 and 150 km/s. The 
angular resolut ion of the map is 9.5". VL A 6cm radio contin- 
uum (contours, lYusef-Zadeh & Morris 1 1987b shows ionized gas 
in the "mini-spiral" near SgrA*. The map is centered on SgrA* 
(red star), (EQ J2000) 17°45'39.9", -29 /, 00'"28.1 I . The two cir- 
cles denote the southern (CND-S) and northern (CND-N) emis- 
sion maxima toward which our excitation studies are performed, 
at offsets of (-20",-30") and (25",40"). The diameter of the cir- 
cles corresponds to the beam of GREAT at the frequency of the 
CO(l 1-10) transition (22.5"). 

2. Observations 

The GREAT observations were performed during two SOFIA 
Basic Science flights in 201 1 . On April 08, GRE AT was operated 
at 427 00 ft altitude in its Lla/L2 configuration dHevminck et al. I 
120121) , tuned to the CO(l 1-10) and CO(16-15) lines. In the sec- 
ond flight, on July 16, (13) CO(13-12), and CO(16-15) were ob- 
served in the Llb/L2 channels (40000 ft). We used the Fast 
Fourier Transform spectrometers (FFTs, iKlein et al71l2012l) for 
back-ends; each provided 212 kHz spectral resolution across 
a bandwidth of 1.5 GHz. At the high observing frequencies, 
GREAT'S bandwidth is insufficient to cover all velocities in the 
Galactic center in one observation. A single setup does cover 
velocities from the individual positions, however, and we made 
two separate receiver tunings offset by +75 and -75 km/s relative 
to rest to cover each position. 

The observations were performed in double-beam switching 
mode, with a throw of 360" at a position angle of 45° (CCW 
to R.A.), at a rate of 1 Hz. We made 5x3 and 3x3 raster maps 
on an 8" grid in equatorial coordinates for CND-S and CND-N, 
respectively. Optical guide cameras established pointing to an 
accuracy of 5". The beam width ©,„/, is 22.5" at the frequency 
of CO(ll-lO), 19.1" forCO(13-12) and 15.6" for CO(16-15). 
The main beam and forward efficiencies are rj m b =0.54 (LI) 
and 0.51 (L2), with rjf =0.95. The data were c alibrated with 
the KOSMA/GREAT calibrator dGuan et al. 1201 2|) , carefully re- 
moving residual telluric lines, and further processed with the 
GILDAS3 packages CLASS and GREG. 

To confine the CO excitation of the CND we have mapped 
during recent years all CO mm/submm transitions, 2 < J up <7 



Transition 



v[THz] 



CND-S 



[Kkm/s] 
CND-N 



12 CO(2-l) 
13 CO(2-l) 
12 CO(3-2) 
12 CO(4-3) 
12 CO(6-5) 
13 CO(6-5) 
12 CO(7-6) 



5.5 

5.3 

16.6 

33.2 

83.0 

79.3 

116.2 



0.230 
0.220 
0.346 
0.461 
0.691 
0.661 
0.807 



797.4 
49.0 
1668.3 
2199.6 
1782.3 
165.5 
1753.0 



670.4 
21.2 
852.9 
1023.5 
786.7 
67.2 
708.6 



I 
1 

A 
A 
A 
A 
A 



12 CO(10-9) 


248.9 


1.152 


841.5 


242.1 


H 


12 CO(11-10) 


304.2 


1.267 


818.8 


218.0 


G 


12 CO(13-12) 


431.3 


1.497 


391.6 


77.0 


G 


13 CO(13-12) 


412.3 


1.431 


26.5 




G 


12 CO(16-15) 


663.4 


1.841 


188.2 


34.4 


G 



http://www.iram.fr/IRAMFR/GILDAS 



Notes. The integration was performed in the velocity ranges [-150,-30] 
and [80,150] for CND-S and CND-N, respectively. All data were con- 
volved to the 22.5" beam of the GREAT CO(ll-lO) observations. For 
each transition we also give the rest frequency and energy of the up- 
per level above ground is. Data are from observations with APEX, the 
IRAM-30m telescope, KWI/Herschel and GREAT/SOFIA. 

with APEX and the IRAM-30m telescope,that are accessible 
from these large ground-based facilities. These data will be pre- 
sented elsewhere (Requena-Torres et al., in prep.); for this study 
we extracted the spectra necessary for a comparison with the new 
SOFIA/GREAT FIR observations. We also included a CO(10- 
9) data point from the HexGal key program on the Galactic 
center (PI: R.Giisten). Table [TJ gives transition information and 
velocity-integrated line intensities for the two positions in the 
CND we analyze here, CND-S and CND-N (see Fig. 1). All 
spectra were convolved from our Nyquist-sampled maps to the 
same angular resolution (22.5", of the CO(l 1-10) observations), 
thereby defining a uniquely coherent data base. 



3. Results and analysis 

Fig. 2 compares the high-J 12 CO spectra observed with GREAT 
with those of lower rotational transitions. At negative veloci- 
ties all low-rotational lines are strongly affected by deep line- 
of-sight absorption features; these are negligible for J up > 6 (E up 
>1 16 K). The FIR lines (J up ^ 10) show clean emission profiles 
from the CND only. To recover the flux affected by these ab- 
sorptions, we fitted a high-order polynomial to the broad emis- 
sion masking the absorptions, and used the fit to estimate the 
underlying emission. Similarly, line contributions from the mas- 
sive GC clouds (at and below +50 km/s) disappear for J up ^10. 
Comparison of the observed isotopic line ratios (~10-20, Table 
1) with the isotopic abundanc e ratio determined for Galactic cen- 
ter clouds, \ l2 C]/\ l3 C] ~2 5 (iGusten et al.ll 1985t IWilson & Rood! 
1994; Ma rtin et aDl2012l) indicates that the CO emission from 
the CND is only marginally optically thick, with opacities on 
the order of unity. 

The line profiles of the high-J CO lines show higher intensi- 
ties for the most negative velocities in the line wings of CND-S 
as the J increases. The rotation curve of the CND shows that 
velocity rises toward SgrA* until just beyond the CND-S and 
CND-N positions, so the increased velocity indicates that the 
higher-J lines are displaced toward the sharp inner edge of the 
CND. While this could indicate heating by the central sources, 
the lack of a strong shift at CND-N suggests that local heating 
is also important. Combined with the broader lines toward the 
brighter CND-S position, this suggests that shocks from turbu- 
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Fig. 2. CO spectra toward the southern (left) and northern (right) 
lobe of the CND. We display the relevant velocity ranges only. 
The data were box-smoothed to 3 km/s spectral resolution. All 
spectra are convolved to the same angular resolution of 22.5". 

lence within the CND m ight provide substantial local heating; 
Christoph er et al.l (|2005) also reported concentrations of dense 
HCN gas cores associated with 2.2 um emission from shock - 
excited H2 at these positions (see, lYusef-Zadeh et aini200ll) . 
Shocks will allow material to drop into the central cavity, while 
radiative excitation will drive material out of the cavity. 

Velocity-integrated brightness temperatures, corrected for 
line-of-sight absorption in multi-component Gaussian fits, are 
compiled in Table 1 . The integration is limited to velocity ranges 
characteristic for the southern [-150 to -30 km/s] and northern 
[80 to 150 km/s] lobe, respectively, thereby excluding velocities 
with heavy line-of-sight absorption and contributions from the 
massive GC clouds. In Fig. 2 the corresponding CO line spec- 
tral energy distributions (SED) are displayed. The peaks of the 
SEDs at J > 7, 155 K above ground reflect the existence of gas 
at high temperatures, and the width of the distribution requires 
contributions from more than a single gas phase. Assuming lo- 
cal thermal equilibrium (LTE), we can estimate the optical depth 
(i"6-5 ~ 1.8, Table 1) and the excitation temperature of the 
bulk emission from th e 12 CO and the 13 CO(6 -5) lines . Toward 
CND-S (clump "O" oflMontero-Castaflo et a l. (2009) and "O" of 
IChristopher et al.l (120051) , with size of ~8") the main beam tem- 
perature T m t, ~50 K observed in the original 8.9" APEX beam 
translates into an R-J corrected clump excitation temperature of 
~135 K. This provides a lower limit for the kinetic temperature 
in this clump. 

The CO excitation was investigated using a large veloc- 
ity gradient (LVG) model in a similar manner as described in 
Wei fi et al71 (|2007). The CO abundances in the GC region are 
well c onstrained by [ 12 CO]/[H 2 ] = 8 • 10" 5 (see iBradford et all 
2005, for details) and, as commented before, the carbon iso- 
topic ratio [ 12 CO]/[ 13 CO] = 25. Furthermore, the velocity gra- 
dient (dv/dr) in the inner CND can be estimated from clump 




1 2 3 4 5 6 7 8 9 10111213141516 1 2 3 4 5 6 7 8 9 10111213141516 

Rotational Quantum Number J upper Rotational Quantum Number J upper 

600 i , , , i^-. ..I m il , . . . , 600 





og(n(H 2 )) [cm 



Fig. 3. LVG model fits with the parameters given in Table 2 for 
CND-S in the left column, for CND-N in the right column. The 
lowest panels display the \ 2 distribution of allowed solutions in 
the density-temperature space (units are from 4 to 14 by 2). 

sizes and line widths observed in high spatial resolution HCN 
and CS interferometry, which yields dv/dr - 150 kms " 1 pc" 1 
(IChristopher et al.l 120051: iMontero-Castano et al.l 120091) . Then 
our model only has three free parameters for each gas compo- 
nent, namely the H2 density, the kinetic temperature T^,,,, and the 
effective source solid angle (filling factor) expressed as equiv- 
alence clump radius ro. In this reduced parameter space the 
line opacities, which are most critical to fit the observed 12 CO 
to 13 CO line ratios, are directly proportional to the H2 density 
(t cc /V(CO/dv) °c[CO]/[H 2 ] n(H 2 ) (dv/dr)- 1 ). As such our mod- 
els have high constraining power on the gas density while T^in 
and ro are partly degenerated. 

The best-fit LVG model parameters to the SEDs are pre- 
sented in Table 2 and displayed in Fig. 2. The main conclusion 
toward the two selected bright rim positions are that 

- the broad, structured SEDs cannot be fitted by a single gas 
component. We must describe the distribution by at least a 
two-component fit as an approximation to a likely continu- 
ous temperature/density distribution. 

In our two component approach the high-J (>10) transitions 
are entirely driven by the high excitation component. This 
component is more pronounced in the south than in the north 
position. Our models suggest that the difference is mainly 
due to a larger source size of this component in the south- 
ern lobe (Table 2). The physical parameters of the two gas 
phases are similar between CND-N and CND-S. The 12 CO 
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line temperatures (with increasing line brightness up to J =4- 
3, Table 1) reflect the lines' low optical depth, driven by the 
high turbulence. The LVG optical depths in the 12 CO(6-5) 
lines are consistent with our LTE estimates above and pre- 
dict optically thin (r < 0.1) emission for the two lowest CO 
transitions. 

- Densities are tightly constrained in the reduced parameter 
space. We find H2 gas densities of log(n) = 4.5*9? cm 4 for 
the low-excitation component for both regions. The gas den- 
sity of the high-excitation component in the southern lobe 
is log(n) = 5.2+9;!. Owing to the absence of a 13 CO(13-12) 
measurement toward CND-N, the parameters for its high- 
excitation gas are less well constrained. Values in Table 2 as- 
sume the same (low) opacity for this line as derived for CND- 
S. Higher gas densities are excluded because they would 
overestimate the line opacities. As such our models strongly 
suggest gravitational unbound clouds. 

- Molecular gas masses (for the restricted velocity ranges, 
Table 1) in the 22.5" beam toward CND-S are M: 3 ° = 280 
(105) M for the lower (higher) excitation phase, and M^p 
= 165 (40) M toward CND-N. Assuming that the CO(7- 
6) emission at all velocities [-150,150] traces CND gas only 
(hence no contributions from 1-o-s gas) the total gas masses 
toward these two positions are 795 and 590 M Q . Next, as- 
suming comparable excitation conditions throughout the in- 
ner CND (defined by an outer disk radius of 3.3 pc), we 
can estimate its total gas mass by scaling the masses de- 
rived in the two beams (1385 M ) with the ratio of the to- 
tal (spatially) integrated flux of the CO 7-6 transition across 
the CND (Fig. 1) to the CO(7-6) flux observed in these two 
beams. We derive a mass of the inner CND of 1.2 10 4 M - 
this is very comparable to the m ass obtained b y Genz eTet al.l 
(1985) and, more recently, by lEtxaluze et al. I (1201 ll) . but 
much lower than results based on dense gas tracers (see dis- 
cussion). 

- Kinetic temperatures for the low-excitation gas are on the 
order of 200 K, and 300-600K for the high-excitation com- 
ponent. The lower limit for Tt,„ of 150 K in all compo- 
nents arises because that lower temperatures would require 
gas densities in excess of the limits discussed above. Kinetic 
gas temperatures as low as the dust temperat ures of 24/45 K 
derive d from recent Herschel observations (Etxaluze et al. 
l201lh are inconsistent with our isotopic line ratios even if we 
consider models with a lower [ 12 CO]/[ 13 CO] abundance ra- 
tio of only 40. Our models therefore strongly support shock 
heating as the most important mechanism because shocks are 
the only mechanism to drive the kinetic temperature signifi- 
cantly above the dust temperature. 



4. Discussion 

The gas densities derived from our CO excitation analysis 
clearly imply that the CO emitting gas cannot be gravitation- 
ally stable against tidal shear because they are much below the 
Roche limit of 10 7 cm" 3 . This comparison and because the ob- 
served velocity gradient is much higher than motions expected 
for virialized clouds at our densities strongly suggest that the 
CND clouds are transient features. This is in contrast to con- 
clusions gained from studies of dense gas tracers such as HCN 
and H CO + (e.g. IChristopher et al.l2005tlMontero-Castafio et all 
2009). These studies derive much higher gas densities (>10 7 
cm -3 ), albeit under the assumption of virialization. As a con- 
sequence, these studies yield molecular gas masses on the order 



Table 2. Parameters of best-fit LVG models 



gas phase r (pc) 



log n(H 2 ) log N(H 2 ) 



CND-S 


low exc. 


0.31 


2001™ 


4 5 +0 - 2 




high exc. 


0.08 


500 +10 ° 


5 2 +0 - 4 
J - -0.2 


CND-N 


low exc. 


0.32 


175 +425 
1 /J -45 


4 5 +03 




high exc. 


0.06 


325 +275 


^ -3+0.6 
-'•-'-0.3 



22.64 
23.34 

22.35 
23.15 



Notes. For each of the two gas phases we quote the kinetic tem- 
perature and H2 density, the equivalent radius (source solid angle) 
and the H 2 column density. The latter is calculated via N(H 2 ) = 
3.09 10 18 n(H 2 ) dv (<5v/<5r) _1 . Recall that model fits are performed in re- 
stricted velocity windows only, Table 1. Uncertainties describe the 3 dB 
X 2 fall-off, see Fig. 2. 

of 10 6 M for the CND, roughly two orders of magnitude above 
our estimate. We note that these differences cannot be explained 
by the different critical densities between CO and HCN/HCO + 
because our high-J CO transitions have values comparable to 
those of the low -J HCN line and even exceed those of HCO + . 

To resolve the obvious discrepancy between these classical 
density tracers, a more rigorous investigation of the HCN/HCO + 
excitation will be performed - including those higher excitation 
submm and FIR transitions that have not been considered so 
far. Complementary 13 C isotopologues shall be observed to con- 
strain the opacity of the lines. Infrared pumping via the strong 
14yum radiation fi eld in the CND will be studied In the ex cita- 
tion models (see ISakamoto et~aTll2010t IWeifi et al. Il2007l for 
an evaluation of the effect). Observationally one should aim to 
constrain this excitation path, and hence the value of HCN to 
trace dense gas in the CND, by searching for the 02 vibrationally 
excited rotational transitions. 

Based on our GREAT/SOFIA data we conclude that the (inner) 
CND is best described by a collection of transient filamentary 
streamers and clumps (Giiste n et al.lll987l) . Its mass is compar- 
atively low, few 10 4 M , which has i mplications on the mass 
accretion rate toward the central object (Gen zel et al.l i2010). 
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